Theory of giant magnetoresistance at misfit interfaces 
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We study theoretically the resistance at the interface between the two planar systems with different 
lattice constants a and b. The resistance and the effect of the magnetic field depends sensitively on 
the ratio a/h. The size of the enlarged unit cell A = uao, = nsfe (nA,ns: integers) is the crucial 
quantity, and the magnetic flux penetrating this enlarged unit cell determines the oscillation of the 
resistance. Therefore, the magnetoresistance is very much enhanced at (nearly) incommensurate 
relation between a and b. 
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Interfaces between different materials are the sources 
of rich physics and functions. Novel phenomena emerge 
that are not expected from each of the constituents [l| . 
One example is a two dimensional metallic state appear- 
ing at the interface between two insulators LaAlOa and 
SrTiOa [ij. Even superconductivity appears in this two- 
dimensional system, which is an issue of recent intensive 
interests Another example is the tunneling magne- 
toresistance (TMR) 0. The resistance across the inter- 
face between the two ferromagnets depends strongly on 
the relative direction of the magnetizations. Therefore, 
transport properties perpendicular to the interface offer 
many useful functions for applications. 

An essential nature of interfaces between different sys- 
tems is the misfit of the lattice constants, which often 
causes the distortion of the lattice structure to relax this 
misfit when it is not so large. In this case, the lattice 
constants slowly change from the interface to the bulk re- 
gion, and correspondingly electrons adiabatically follow 
this gradual change. When the misfit is larger, on the 
other hand, the system can not remedy this misfit and 
an incommensurate situation occurs at the interface. 

Incommensurate systems attract recent attention from 
the viewpoints of charge/spin density waves 4|, localiza- 
tion of wavefunctions [5| , and quasi-crystals 6 1 . In these 
systems, incommensurability occurs in the bulk states, 
which is rather exceptional or special cases. On the other 
hand, the incommensurability occurs very often at inter- 
faces since there is no definite relation between the lattice 
constants of the two systems. 

In this paper, we study theoretically the tunneling 
conductance across the interface between the two two- 
dimensional systems A and B with lattice constants a 
and 6, respectively. The ratio a/h matters significantly 
in this tunneling process, and also its sensitivity to an 
external magnetic field. 

Let us start with the two one-dimensional chains A and 
B. The extension to two-dimensions is straightforward. 
We assume two chains have the same length L. The 
number of sites in chain A and in chain B are Na and 
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FIG. 1: The upper panel indicates two chains in the real 
space. One chain has the lattice constant a, and the other 
chain has the lattice constant b, which are placed in a mag- 
netic field. The smallest periodic part is constructed from 
riA sites of chain A and ub sites of chain B, i.e., the size of 
the enlarged unit cell A is given by A = an a = bnB- Gab is 
the number of this enlarged unit cell and hence is the great- 
est common divisor of Na and Nb, i.e., Na = nAGAB and 
Nb = nBGAB, and L = GabX. The bottom panel indicates 
the momentum space of chain A and chain B. The Brillouin 
zones of A and B are discretized by the same unit 2-k jh, while 
the sizes are different. The Brillouin zones are decomposed 
into riA and ns parts by 



Nb-, which are determined from L — uNa — bNs- We 
assume that two chains are parallel to z-direction with 
the separation c, and the tunneling amplitude between 
the nth site in chain A and the mth site in chain B is 
given by 
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(1) 



Here d characterizes the spatial extent of the tunnel- 
ing process. Two chains are placed in a magnetic field, 
which is perpendicular to the plane including two chains. 
The magnetic field induces AB phase (Aharanov-Bohm 
phase). We choose the gauge as A = (0, Bx). 
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We rewrite Eq. ([l} by wavenumber representation as 
1 



with 



(2) 



where the wavenumbers are specified by the integers k 

and p. The lattice constant A of the composite system of 

A and B is given by A = an a — bns, where we define as 

71^1= J^-^ and ub — ■ Gar is the greatest common 
, , (jab (jab ^ ^ 

divisor of Na and Nb^ and is the number of the unit 
cells with the lattice constant A, i.e., L = Gab^- The 
translational symmetry by Xe^ leads to the conservation 



of wave numbers by mod 



The Brillouin zones are 



decomposed into ua and tib parts by 

The summations in Eq. ([2]) can be carried out (See in 
the appendix) and tkp is obtained as 
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cosh i - cos teifc-^^, 



6gab is defined as 



Sgab {k-p) 



1 k — p = mod Gab 
otherwise. 



(5) 



Equation ([5]) represents the wavenumber conservation. ^ 
is the characteristic length of this system, which is repre- 
sented as £ = — — _2__ xa and xb are solutions 

^ ^ Lab ua riB -° 

of the Diophantine equation UbXa — nAXB = 1- This 
equation has an integer solution because ua and ub are 
coprime. (xa + mnA,XB + muB) is also the solution, 
where m is an arbitrary integer, but this indefiniteness 
is not concerned with Eq. ([3]) because of Eq. ([S|) . Equa- 
tion (|4]) represents the interference in the periodic part. 

We consider the resistance at an interface between 
the two-dimensional systems A and B, which are the 
straightforward generalization of the above chain system. 
The conductance g per one enlarged unit cell is repre- 
sented as 
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We take a limit Gab oo, and 
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Fik,p) = \tABfc{k..p.)fciky,Py)r (8) 
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b is the dimensionless magnetic flux penetrating the en- 
larged unit cell, i.e., BcX — 27r6. The resistance R per 

,2 

the unit area is given by i? = which is the physical 
quantity of our main interest. 

First, we study the lattice constant A of the enlarged 
unit cell, j change radically as the ratio f = ^ — 
Particularly, j indicates a fractal architecture with 
flxed Na, if iV^ is a power of a prime number. We 
show this relation in Fig. [2] with the flxed lattice con- 
stant a of chain A. The upper left panel indicates the 
global behavior of j v.s. | = and the rest pan- 
els show the graphs for selected ^ with fixed Na — 
2l",3^64. When Na = 2^0, 3«, it shows clear frac- 
tal structures, while it shows a complex structure when 
Na — 6*- When Na — p" where p is a prime num- 
ber, J is given by the following calculations. We define 
Si = {p"-%2p"-%---(p*-l)p"-'} {i = l,2,---n) and 
Pi — Si — Si-i [i — 1,2, •••«), where we set Sq — (j) 
(empty set). The relation between A and is specified 
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when Nb G Pi- This relation and the definition of Pi 
generate the fractal structure scaled by i. When Na is 
not a power of a prime number, this fractal structure is 
not there, but shows a highly singular behavior as a 
function of | as shown in panels (a) and (d) of Fig. [2] 

Next, we consider the magnetic field dependence of 
R to see the proper scaling for B. In a limiting case 
d — >■ 0, the hopping amplitude is finite only between sites 
whose x-coordinates and y-coordinates are the same. In 
this limit, F = 1 and it is clear that the resistivity has 
the period /Shi = 1, which is independent of {nA,nB)- 
From BcX = 27r6, the resistivity of the misfit interface 
is enhanced for larger A (or ua, ub), i.e., (nearly) incom- 
mensurate case. In this case, f = — is small, so we 
cannot assume ^ 3> c?. Then there is a crossover from 
(i— >0<cOof^^fi and the interference in the peri- 
odic part F occurs in the latter case where the period of 
the conductivity becomes A6i = 21 ab, where Iab is the 
least common multiple of ua and ub- However, as will 
be seen for explicit examples, the variation of R occurs 
within the scale of A5i = 1 even for c? ^ ^. 

We estimate the behavior of R with A in the limiting 
cases. In a limiting case d — ^ 0, one enlarged unit cell 
has one pair of sites connected by the finite hopping am- 
plitude, and hence R grow as cx A^. When d = a = b, 
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FIG. 2; The variation of j (A the size of the enlarged unit 
cell) as a function of ^ = ^ : Panel (a) indicates the global 
behavior while the rest panels (b),(c), and (d) indicate j for 
selected f = ^ with fixed Na = 2^°, 2,^, and 6*, respectively. 
When Na ~ 2^", 3®, it show fractal structures, but it shows a 



complex structure when Na = 6 



which is more relevant to the reahstic systems, there are 
many finite hopping amplitudes between sites in one en- 
larged unit cell, and the number of site in each chain does 
not becomes important. The hopping amplitude per area 
does not change with A and the scale of R is nearly con- 
stant although the magnetic field dependence is sensitive 
to A. In the extreme case a <C 6, on the other hand, the 
situation is similar to the case of d — 0, and R is expected 
to grow with oc A^. 

Now, let us study a concrete example. We assume the 
simplest tight binding model on a square lattice both for 
A and B as ^a(b) (k) — cosk^ + cos ky with the chemical 
potential at /i = 0, and the magnetic field is along the 
X-axis, i.e., by = 0. From Eq. (|6]) and Eq. ([7]), we can 
obtain the analytic form of the conductivity as given in 
the appendix. In the following calculations, we fix the 
lattice constant of ^ as a = 1 and set = 1. 

In this model, the Fermi surfaces is straight line as 
indicated in Fig. [31 At the misfit interface, the Brillouin 



zone is reduced to nj^ x {ng x ub) parts by The 
segmented pieces of the Fermi surfaces are determined by 
the parity of n^i (n^). This even-odd effect is a special 
properties of the linear Fermi surfaces. In Fig. [3l we 
indicate two examples, i.e., cases of ha — 2 and 3. 

We now consider the resistance i? as a function of the 
magnetic field iJ^ and the dimensionless magnetic fiux 
bx- The resistance R a.t = b^ = reflect the even-odd 
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FIG. 3: The Fermi surface for €a{b) (fc) ~ cos fcj, +cos fcj,: The 
Fermi surfaces are indicated by red lines. Due to the folding of 
the 1st Brillouine zone (BZ), the Fermi surface is segmented 
into smaller pieces in the reduced BZ, which depends on the 
parity of ua (wb). We show two examples ua ~ 2 and 3 at 
the right panel for even and odd nA(ns), respectively. 
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FIG. 4: Resistance per unit area R for d — 0: The panels (al) 
and (bl) shows i? as a function of bx while (a2) and (b2) as 
a function of ^f^- Panels (al) and (a2) are for (nA,ns) = 
(2,1), while (bl) and (b2) are for (nA,ns) = (50,49). R 
depends on the magnetic field through the flux penetrating 
the enlarged unit cell, and when A is large, R oscillates as 
more rapidly and the maximum value becomes larger as 
oc A^ as shown in Fig.6(bl). 



effect of the Fermi surfaces. When both ua and ub are 
odd, i? = for integer values of bx- On the other hand, 
when either ua or ub is even, R — at bx — integer -I- ^. 
Strictly speaking, R — means the diverging g and hence 
the perturbation theory with respect to Iab does not 
work there. This divergence has two origins, i.e., the 
Van Hove singularity and the fact that the forms of the 
Fermi surface of eA and es are the same. 

At first, we consider a limiting case — > 0. The be- 
havior of the resistance R depends on the parity of ua 
and Ub, which reflects the even-odd effect of the Fermi 
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FIG. 5: Resistance per unit area R for d — a: The panels 
(al) and (bl) shows _R as a function of bx while (a2) and (b2) 
as a function of Panels (al) and (a2) are for {nA^ns) = 

(2, 1), while (bl) and (b2) are for {nA,nB) = (50,49). 
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FIG. 6: The maximum value of the resistance R in the range 
of < B^ac < 27r as a function of f = ^ for d = ((al)) 
and d — a ((a2)), respectively. This complex relations are 
cleanly organized as a function of A, which are shown in (bl) 
and (b2). The black solid lines are the guide to the eyes (slope 
2) for the asymptotic relation R (x \^ . 



surfaces. The maximum value of R changes radically in 
each nA and ub, which will be discussed later. R has 
the period from bx = to = 1, and indicates a similar 
behavior to bx for all (n^,ns) except for the even-odd 
effect. When either ua oi ub is even, the bx is shifted 
by ^. bx is the dimensionless magnetic flux penetrating 
the enlarged unit cell, which is scaled in an appropriate 
manner in each case as BxcX = 27rbx- The larger A be- 
comes, R oscillates more rapidly as Bx- The panels (al) 
and (bl) of Fig.|4]show i? as a function of bx, while (a2) 
and (b2) shows i? as a function of Panels (al) and 

(a2) are for {ua, ub) = (2, 1), while (bl) and (b2) are for 
[uatUb) = (50,49). If A become bigger, the resistivity 
R oscillate as more rapidly. 

Next, we consider a finite d = a. The behavior of the 
resistance R is not determined only by the parity of ua 
and Ub, but it is different in each ua and Ub because of 
the interference term F. The periodicity of R changes 
from M)x = 1 to ISbx — 21 ab, but the 6a;-value at which 
R = remains unchanged from the case of d = 0. The 
panels (al) and (bl) of Fig. [5]show i? as a function of bx, 
while (a2) and (b2) as a function of Panels (al) 

and (a2) are for {nA,nB) = (2,1), while (bl) and (b2) 
are for [uatTLb) — (50,49). In the similar way to the 
case oi d = 0, the larger A become, the R oscillates as 
bx changes of the order of 1 and hence more rapidly as a 
function of 

Establishing the enhanced magnetoresistance by A, we 
next discuss the maximum value of the resistance i^max 
in the region < Bxac < 27r. -Rmax indicates a singular 



behavior as a function of | = ^ in (al) and (a2) of 
Fig. |6] for two cases d = and d — a, respectively. Both 
show the rather complex and singular behavior, but these 
are neatly organized as a function of A as shown in (bl) 
and (b2) of Fig. (6] When d -> (((bl)), i?max increases 
as A as expected for all the regions of a/b (asymptotically 
A^ in the large A limit). For d — a (((b2)), on the other 
hand, i?max stays almost constant and independent of A 
for a = b, while it approaches to the behavior of d = as 
a/b decreases. 

Now we discuss about the relevance of the present re- 
sults to real systems. The disorder effect at the inter- 
face gives the mean free path £. When the size A of the 
enlarged unit cell is larger than £, the singular depen- 
dence on a/b is broadened. In other words, the enhance- 
ment of the magnetoresistance saturates by the factor 
^min(^, A)^. The most relevant case to the real systems 
is that a = b and d = a. In this case, the scale of the re- 
sistance R per unit area does not sensitively depends on 
the ratio a/b, while the magnetoresistance is determined 
by A and depends strongly on the ratio a/b in a singu- 
lar way. The essence of the enhanced magnetoresistance 
is the sensitive change in the interference pattern of the 
wavefunctions within the enlarged unit cell induced by 
the magnetic flux, it is expected that the magnetic field 
perpendicular to the interface also gives the similar effect 
to the parallel case discussed in the present paper. 

In summary, we have studied the magentoresistance at 
the interface with misfit of lattice constants. We found 
that resistance R depends on the ratio a/b of the two 
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lattice constants in a singular way, and the size A of the 
enlarged unit cell determines the magnitude of the mag- 
netoresistnce, which can he enhanced orders of magni- 
tudes when a/h is a (nearly) irrational number. 
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detail calculations of Eq. (3) 



In this appendix, we show detail calculations of Eq. (3). We write Eq. (1) by wave-number representation as 

U ^ e"'^"+*^'" ill) 



nm 



tAB 



nm 



, where we use a = i^qj^ and b = :^gj^- first, we decompose the summations into ha and hb parts by Gab, 

^/NaNb „^„^=0 3,^0 y^Q 

I Gab— 1 nA— InB — 1 % / % 

JNaNr ^ ■ \ } 



We define j/_ = maxjy : UbX > n^y} and rearrange the order of the summations, which begins at y_, by making use 
of periodicity. It is specified as 

tAB ^^-^ ,-#-X^(»BX-..fa+l+.-))-. -("'^7+-) (fe-M)+. -("'-X-"^""^ f 16) 

V-^Ai-^B „^=o x=0 j/=0 

Here, the summations over ruA and y is easily carried out. tkp becomes 

^to = ^AB^^"" (fc-p-2M) 1-er^ ^^(„,(,_+i)_„,.)_,^(,_M)+iM^(p+M) , . 

t, J''-^^^ (/v-i>-2M) l-<-^ ^_^(^r^B<c-nA{l+V-))-i^{k-M)+ ^<l-^+'\ p+M) 

We define I as I = ~ T^Z ~ V- '^^ equal to the quotient of ubx divided by ua- We denote the remainder as 
Aa;, i.e., ubx = uaV- + A^^. We write tkp by Aa; and k — p — 2M = uGab, 

_ 6GABik~p-2M) 1-e-? ^^+,^(p+M)"v^' -»|ff»+(i-i^(p+M))Ax , . 

_ jGABitzl::^ 1-^"^ el+^llt(P+M) -.a^n-(l+.^(p+M))Ax , . 

We define xa and a:s as solutions of the Diophantine equation, which is expressed as 

ubxa - uaxb = 1- (21) 
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This equation has a integer solution because ua and ub are coprime. The general solutions of this equation is 
represented as {xa + mnA,XB + "muB), where m is an arbitrary integer. Thus Xa = Axxa + muA and Xb = 
Axxb + muB are satisfied with 



ubXa - uaXb = Ax. (22) 

There is one-to-one correspondence between x — 0, 1 • • -ua — 1 and Ax — 0, 1 • • -ua — 1 because ua and ub are 
coprime. Ax is represented as Ax = Axxa + muA where m is dependent on x. The summation over x is transformed 
into the summation over Ax, 



After the summation over Ax and some calculations, tkp become 



^ ^ cosh (1) - COS (^(p + Af) + ffx^) 

Eq. (3) is given by arranging Eq. 1251 into the symmetric form by using k ~ p — 2M = hGab, which is 

= -4^Sa.^ {k~p- 2M) (l-e-i) ^Ml) (26) 



detail calculations of Eq. (6) and Eq. (7) 



In this appendix, we show detail calculations of Eq. (6) and Eq. (7). We assume that the tunneling amplitude 
between two lattices is the multiple of two copies of Eq. (1). The wave-number representation of the tunneling 
amplitude also become the multiple of two copies of Eq. (3), because calculations can be carried out independently 
in X- and y-directions. It is specified as 



f2D {k-M,p + M) ^2 ^ 

ikp^^AB cIgab ("^ ^ P ^ 2M) (27) 

UATlB 

f2Dik,p) = f{k,,p,)f{ky,py) (28) 

^Gab (fe - P) = ^Gab (kx - Px) Sgab {ky ^ Py) ■ (29) 



The conductance between two lattices is calculated from 



G^j:h,fs(.4^^t))sL(^^p)). (30) 



kp 



^A(B) is the energy dispersion for lattice A(B), which are periodic by 27r. The Brillouin Zones are similarly decomposed 
into Ua X riA parts and ub ^ ns parts by new wave- number conservations by translational symmetries by Xe^ and 
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Agj,. Eq. (j30p is transformed into the more meaningful form by the fohowing calculations. 

2 Na-1 Nb-1 / /r, w / /r, 



-4 B k^,kv=Op^,Py=0 
''AB 



E E l/2D(P + M + nG^B,P + M)|'^(^^(^(nGAB+p + 2M)jj5Ui3f^P 



,2 "A — 1 Gab — 1 ns — l 

''AB 



E E E \f2D{GABm + i + M + GABn,GABm + i + M)y 



9 2 

,ny — ,2y — -''■'T-y — 

X S (^U (J^ {GABn + Gabtu + i + 2M)j j S (^£,b {GABm + i) 

2 "A- l G ab- 1 »b- 1 

E E E l/2i3(GABn + i + M,GABm + i + M)|' 

-'^ ,ny — ,2y — rria: ,my — 

E E <^«(M) 

^4 71^ ,ny— ma; -TMy— 



with 

,2 Gab-1 



= E \f2D{GABk + i + M,GABP + i + M)\'5{^A{GABk + i + 2M))6{iB{GABP + i)) (31) 

Here, we take a limit Gas oo. Each term is replaced by -^^i 2ttx and -4^Af — > 27rb, 



^nAnB(^)^ ll^ j'xF (^k 



Gab 

X + \,p + X+'^5UA{^^{k + b + x)\\8UB (^iP + -) 



with 



F{k,p) = \tABfc {k^,Px) fc iky^Py)f 
sinh (i,) 



cosh (i)- cos (^fc-^p) 
/2 = [0,1] X [0,1]. 

detail calculations of the concrete example 



In this section, we calculate the conductivity of the concrete example, i.e., £a{b) (k) — cosk^ + cosfcj, and by — 0. 
The conductivity is calculated from 



E E /L^=^^ .A — (fc + 6e. + =r) 5 eB — + 



KA — 1 "B-l 



j2 '"'A'^B \ \nA / / \ Kf^B 



The calculation is done by the simple and straightforward method. We carry on the delta functions by the rule 
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one by one, where x, is the ith zero point of g {x). The result is 

Ann — ~. o 



kx=0 Px=0 







-"^ ' 2 ' 2 ' 2 ' 2 



Sill (fcc. +b. + A)) sin (p, - 6, + A)) 
with A = — J . In the limit d — )• 0, The conductivity is summarized into the simpler form as 



(33) 
(34) 



n^— 1 riB— 1 



\ y y 

"^^^^AnB I sin (fc, + A + 6,)) sin (p, + A - 6,)) 



(35) 
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